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Values for the bulk viscosity r)b of molecular nitrogen gas (N2) were derived from spontaneous Rayleigh- 
Brillouin (RB) scattering at ultraviolet wavelengths (A = 366.8 nm) and at a 90° scattering angle. Analysis of 
the scattering profiles yield values showing a linear increasing trend, ranging from r] b = 0.7 X 10~ 5 to 2.0 X 10~ 5 
kg-m _1 s _1 in the temperature interval from 255 K to 340 K. The present values, pertaining to hypersound 
acoustics at frequencies in the GHz domain, are found to be in agreement with results from acoustic attenuation 
experiments in N2 performed at MHz frequencies. 
© 2013 Optical Society of America 
OCIS codes: 010.0010, 290.5830, 290.5840, 290.5870. 
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The concept of bulk viscosity, 775, also referred to as 
volume viscosity is part of a thermodynamic description 
of gases as a transport coefficient in addition to the shear 
viscosity i] s [1,2]. The bulk viscosity results from the col- 
lisional energy exchange between the translational and 
internal (rotational and vibrational) degrees of freedom 
in fluids. The value of 775 of gases can be measured via 
sound absorption but only a limited number of studies 
have been reported [3,4]. Furthermore, such measure- 
ments yield values for r/i, related to acoustic frequen- 
cies in the MHz range, while bulk viscosity is regarded 
as a frequency-dependent parameter [5], resulting from 
the competition between the internal relaxation time of 
molecules and the period of acoustic waves. Therefore, 
the values measured at MHz frequencies should not be 
directly applicable to much higher frequencies, such as 
in light scattering experiments, where the (hypersound) 
acoustic waves are in the GHz domain. For example, 
Pan et al. found that bulk viscosity for CO2 in their 
coherent Rayleigh-Brillouin scattering (CRBS) experi- 
ment, is 1000 times smaller than the sound absorption 
value [6] . They suggested that values of bulk viscosity at 
high frequencies could be derived by comparing the light 
scattering profiles of gases to accurate models developed 
by Bolcy et al. [7] and Tenti et al. [8] , given that in these 
models the only unknown parameter is 

In this Letter, we present measurements of sponta- 
neous Rayleigh-Brillouin scattering (SRBS) profiles of 
N2 in a temperature range of 255 K to 340 K and a 
pressure range of 850 mbar to 3400 mbar. The measured 
scattering profiles are compared to the so-called Tenti 
S6 model [8], which is generally considered as the most 
accurate model to describe the RB-scattering profile [9] . 
Implicit in the model is that the Brillouin side peaks to 
the central Rayleigh peak in the scattering profile are 
shifted by [10]: 



with n and v the index of refraction and the sound veloc- 
ity in the gas, and ui and 9 the angular frequency of the 
light and the scattering angle. The Brillouin side peaks 
exhibit a profile, associated with the damping of acoustic 
waves, and dependent of the thermodynamic properties 
of the gaseous medium as well as the light scattering pa- 
rameters, yielding a Lorentzian profile of full width half 
maximum: 
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p the density, k the thermal conductivity, and 
7 = Cp/C v . The code implementing the Tenti model 
(version S6) was based on that of Pan et al. [11], and 
was used for previous studies on spontaneous and coher- 
ent RB scattering in gases [12, 13]. This method via the 
Tenti model must be followed for extracting rjb in gases 
where the central Rayleigh peak overlays the Brillouin 
side peak, unlike for liquids where the Brillouin features 
are fully isolated and 77^ can be determined directly by 
measuring the width Tb [14]. 

Details of the experimental setup and methods for 
measuring high signal-to-noise RB scattering profiles 
have been reported in [15]. The profiles are recorded 
for scattering at 9 = 90° induced by an effective intra- 
cavity circulating power of 5 Watt at A = 366.8 nm, via 
a plano-concave Fabry-Perot interferometer (FPI) with 
an instrument linewidth of 232 MHz. For each measure- 
ment, the scattering cell is initially charged to one of the 
designated pressures, namely 1 bar or 3 bar, at room 
temperature, followed by sealing the cell and then set- 
ting the temperature to one of the designated values: 
255 K, 275 K, 297 K or 336 K. The actual pressure of 
each measurement thus differs from the initial pressure, 
while the number density of the gas molecules remains 
the same. The actual pressure is derived via the ideal gas 
law. 

Scattering profiles of N 2 at eight different (p, T) 
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256.0 K 0.862 bar n = 0.75x1 0" 5 275.0 K 0.940 bar q = 1 .05x1 0" 5 296.7 K 1 .01 bar r|„= 1 .4x1 0" 5 336.0 K 1 .140 bar n = 2.0x1 0" 




frequency (GHz) frequency (GHz) frequency (GHz) frequency (GHz) 



Fig. 1. Rayleigh-Brillouin scattering profiles (black dots) as measured for various (p, T) pressure-temperature com- 
binations as specified. A comparison is made with calculations via the Tcnti S6 model (red lines), convolved for the 
instrument width of 232 MHz and for values of the bulk viscosity i]b, deduced from the profiles recorded at 3 bar. 
Residuals between the measurements and the calculations are given underneath. 



pressure-temperature combinations are shown as black 
dots in Fig 1. Since the effect of r\b is most significant 
at the highest pressures, where the Brillouin side peaks 
become pronounced (see Fig. 1), the data recorded with 
an initial pressure of 3 bar are used for determining r\b- 
Fig. 2 elucidates the method for extracting a value for 
7]b in the comparison of the Tcnti S6 model with the RB- 
profiles for the specific recording of an RB-profilc in N2 
under conditions T = 336.6 K and p = 3.40 bar. Panel 
(a) shows the measurement (black dots) and the mod- 
eled scattering profiles, for three different values of bulk 
viscosity, and for values of the N2 transport coefficients 
as obtained from literature (listed in Tabic 1). For the 
dimcnsionlcss internal specific heat capacity of internal 
degrees of freedom Ci„t a value of 1 is used throughout. 
Residuals between the measurement and the three mod- 
eled scattering profiles arc shown in (b). Panel (c) shows 
a ^-calculation as a function of bulk viscosity employed 
in the Tenti S6 model. 

Table 1. Transport coefficients used for modeling the RB- 
profiles of N 2 . Values for rj s and k are calculated accord- 
ing to the Sutherland formula in [16], and rib from the 
present experiment. 
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This procedure of optimizing rjb was applied to the 



RB-scattering measurements for initial pressure of 3 bar 
N2. The resulting values for r/b and their uncertainties 
are plotted in Fig. 3, combined with values from litera- 
ture. Prangsma et al. [4] determined bulk viscosities for 
N2 using sound absorption measurements in the tem- 
perature range T = 70 ~ 300 K. The experiment inves- 
tigated a wide range of acoustic frequency to pressure 
ratios, but all in the MHz domain. Pan et al. [11] used 
the value from Prangsma et al. [4] and found good agree- 
ment between their CRBS profile and a calculation using 
the Tenti model (the S7 variant) [7], suggesting that the 
value of bulk viscosity for N2 obtained at MHz frequen- 
cies is also valid for the GHz range. Cornelia et al. [17] 
successfully modeled CRBS profiles in N2 assuming a 
r)b/r) s ratio of 0.73 from [4], valid at room temperature, 
and extrapolated this to 500 K. Values previously ob- 
tained by Vietez et al. [12] using SRBS at 3 bar N2 
slightly deviate; however, no uncertainty was specified 
and if a similar uncertainty is assumed as in the present 
study agreement within combined la follows. Mcijer et 
al. [18] using CRBS (at 532 nm) at 5 bar N2 deduce an 
even larger value, but still agreement within 2a results. 

The present experimental results for 775 in the tempera- 
ture interval 254 — 337 K, shown as black dots in Fig. 3 
show a linear dependence with temperature, roughly 
similar as in [4]. While the data of [4] extend to tem- 
peratures as low as 180 K, and the present data extend 
to 337 K, for the overlapping range 250 — 300 K good 
agreement is found. It is assumed that for dilute gases 
the bulk viscosity is independent of pressure, similar to 
shear viscosity and thermal conductivity [16]. The RB- 
profiles recorded for 1 bar N2 gas, shown in the upper 
panels of Fig. 1, are modeled with the rjb(T) values ob- 
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Fig. 2. (a) Experimental Raylcigh-Brillouin scattering 
profile in N2 for 3.40 bar and 336.6 K (black dots), and 
convolved Tcnti S6 calculations for bulk viscosity being 
1.0 xf(T 5 (green line), 2.0 xl0~ 5 (red line) and3.0xl0~ 5 
(yellow line) kgm _1 s _1 , respectively, (b) Residuals be- 
tween measured and calculated scattering profiles for 
these three values of the bulk viscosity, (c) A plot of the 
\ 2 as a function of bulk viscosity. The optimized value 
of bulk viscosity is found at the minimum of % 2 , with 
the gray area indicating the estimated statistical error, 
calculated according to procedures discussed in [12, 18]. 



tained for 3 bar, also yielding good agreement. While 
the shear viscosity rj s is known to exhibit a linear tem- 
perature dependence in the window 254 — 337 K [16], 
the ratio r]b/Vs grows from 0.46 to 1.01 for the present 
data. This behaviour may be related to the freezing out 
of internal degrees of freedom at lower temperatures. 

A general conclusion is drawn that for pure nitrogen 
(N2) gas the bulk viscosity at acoustic frequencies in the 
MHz regime is the same as for hypersound frequencies in 
the GHz regime. This result is surprising in view of the 
results in carbon dioxide (CO2) gas where differences by 
orders of magnitude were found [6]. 
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Fig. 3. Comparison of bulk viscosity measured from dif- 
ferent experiments. Note that the result of Pan et al. [11] 
overlays a data point by Prangsma et al. [4]. Data of 
Vieitez et al. [12] and Meijer et al. [18] also included. 
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